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TRANSVERSE  BEAM  DYNAMICS  IN  THE  MODIFIED  BETATRON 


I.  INTRODUCTION 

It  has  been  suggested1'2  that  the  current  carrying  capacity  of  a  conventional  betatron  accelerator 
might  be  improved  dramatically  by  the  addition  of  a  strong  toroidal  magnetic  field.  Such  a  field  acts  to 
confine  the  beam  during  injection  and  early  stages  of  acceleration  when  y ,  the  usual  relativistic  factor, 
is  small  and  space  charge  effects  which  tend  to  expand  the  beam  are  large.  After  acceleration  is  com¬ 
plete,  y  is  large,  space  charge  effects  are  small,  and  the  usual  weak  focussing  betatron  fields  are 
sufficient  to  confine  the  beam;  the  toroidal  field  may  then  be  removed  to  facilitate  beam  ejection.  In 
general  both  vertical  and  toroidal  magnetic  fields  may  be  changing  simultaneously  during  beam  injection 
and  ejection.  It  is  the  purpose  of  this  paper  to  examine  the  behavior  of  the  beam  in  such  time-varying 
fields. 

We  shall  derive  and  solve  equations  governing  the  motion  of  the  center  of  an  electron  beam 
confined  in  a  modified  betatron  as  well  as  equations  governing  the  motion  of  an  individual  particle 
within  the  beam.  Whole  beam  and  single  particle  stability  criteria  will  be  presented;  (he  stabilizing 
effect  of  the  toroidal  field  for  both  beam  and  single  particle  motions,  noted  earlier,1  2  will  be  apparent. 

When  the  fields  are  allowed  to  vary  in  time  two  interesting  phenomena  occur.  The  first 
phenomenon,  which  occurs  during  acceleration,  has  no  analogue  in  a  conventional  betatron:  As  the 
beam  accelerates  (y  increases)  the  betatron  makes  a  transition  from  a  region  in  parameter  space  in 
which  the  toroidal  field  is  essential  to  stability  (modified  betatron  regime)  to  a  region  in  which  the 
toroidal  field  is  superfluous  to  stability  (conventional  betatron  regime).  It  turns  out  that,  except  under 
extraordinary  circumstances,  the  system  must  pass  through  an  "instability  gap”— a  region  of  parameter 
space,  separating  the  modified  and  conventional  betatron  regimes,  in  which  single  particle  motion  is 
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unstable,  though  beam  center  motion  may  not  be,  irrespective  of  the  magnitude  of  the  toroidal  magnetic 
held.  However,  though  the  size  of  the  instability  gap  is  independent  of  the  toroidal  field,  the  instability 
growth  rate  within  the  gap  is  inversely  proportional  to  this  field.  We  find  below  that  by  judicious  mag¬ 
net  design  and  sufficiently  rapid  acceleration,  this  gap  may  be  successfully  traversed  with  minimal  beam 
disturbance. 

The  second  phenomenon  occuring  in  time  varying  fields  does  have  an  analogue  in  a  conventional 
betatron;  this  is  the  adiabatic  change  in  the  amplitude  of  the  betatron  oscillations.3  Since  the  frequency 
of  these  oscillations  depends  now  on  both  the  vertical  and  toroidal  fields  a  slow  change  in  either  is 
expected  to  alter  the  amplitude  of  the  betatron  oscillations.  During  acceleration  we  find,  as  in  a  con¬ 
ventional  accelerator,3  that  the  oscillation  amplitude  decreases  as  the  vertical  field  increases.  If  one  now 
considers  removal  of  the  toroidal  field  prior  to  beam  ejection,  we  find  that,  as  long  as  the  toroidal  mag¬ 
netic  field  is  much  larger  than  the  vertical  field,  the  beam  motion  will  describe  orbits  of  increasing 
amplitude  as  the  toroidal  field  is  decreased.  Once  the  toroidal  field  becomes  comparable  to  the  vertical 
field,  however,  the  motion  becomes  more  complicated  and  the  betatron  oscillations  no  longer  continue 
to  increase  in  amplitude.  We  find  that,  by  careful  choice  of  field  strengths,  the  ratio  of  the  betatron 
oscillation  amplitude  before  acceleration  to  the  amplitude  of  oscillation  following  complete  removal  of 
the  toroidal  field  can  be  adjusted  to  be  near  one. 

In  the  following  analysis  we  assume  "perfect,”  i.e.  azimuthaily  symmetric  fields.  By  neglecting  the 
possibility  of  azimuthal  variation  in  the  self  fields  (due  to  beam  bunching  or  kinking)  we  omit  here 
consideration  of  a  variety  of  beam  instabilities  that  may  occur;4  by  neglecting  similar  azimuthal  varia¬ 
tion  in  the  applied  fields  ("field  errors")  we  neglect  the  effects  of  orbital  resonances.  These  will  be 
addressed  in  a  separate  report.5 
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II.  EQUILIBRIUM  RADIAL  FORCE  BALANCE 

The  geometry  of  the  modified  betatron  is  shown  in  Fig.  1.  The  field  configuration  is  that  of  an 
ordinary  betatron  with  the  addition  of  a  toroidal  magnetic  field,  B„ 0,  here  taken  to  be  positive  and  con¬ 
stant  across  the  minor  cross  section  of  the  torus.  We  consider  an  electron  beam  of  circular  cross  sec¬ 
tion,  as  shown  in  Fig.  2,  with  center  located  at  (rc,  zf)  -  (r0  +  Ar,  Az)  where  r0  is  the  equilibrium 
radius  for  the  center  of  the  beam  at  which  the  electric,  magnetic,  and  centrifugal  forces  on  a  particle  at 
the  center  of  the  beam  are  in  balance.  We  shall  take  r0  to  be  the  major  radius  of  the  accelerator 
chamber.  In  the  absence  of  self  field  effects  radial  force  balance  requires  the  electron  circulation  fre¬ 
quency  at  r  -  r0,  z  -  0  to  be  given  by 

0O  ”  fl  -o  =  eBt0/my0c  (no  self  field  effects)  (1) 

where  B:0  is  the  value  of  the  applied  vertical  betatron  field  at  the  location  of  the  orbit,  y0  is  the  usual 

relativistic  factor,  e(>0)  is  the  magnitude  of  the  electron  charge,  m  is  the  electron  rest  mass,  and  c  is 

the  speed  of  light. 

Self  field  effects  will  modify  Eq.  (1)  however.  A  nonneutral  current  ring  produces  both  a  zero 
order  vertical  magnetic  field  and  a  radial  electric  field.  In  general,  for  a  reference  particle  at  r  -  r0,  z  = 
0,  radial  force  balance  requires 

-yo'>o  "  +  7'oM;,0,|  (2) 

where  Er{0)  and  B.i0)  are  the  zero  order  fields  at  r  -  r0,  z  -  0.  From  Appendix  A,  Eqs.  (A-25c,  26c, 
26d) 

B}m  "  B;q  —  n  n<)efio  lB  (3) 

ro 

Eri0)  -  -ir«oC—  l£  (4) 

ro 

where  the  notation  is  defined  in  Appendix  A. 

The  terms  proportional  to  lg  in  Eq.  (3)  and  l£  in  Eq.  (4)  are  toroidal  corrections  to  the  self  fields 
of  a  cylindrical  beam.  They  represent  "hoop  stresses?  —  self  forces  on  a  nonneutral  ring  of  current 
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which  act  to  expand  the  ring.  Since  we  do  not  attempt  here  to  construct  a  consistent  equilibrium  for 
the  beam  we  leave  lB  and  lE  arbitrary  in  the  analysis  below  since  their  precise  values  depend  upon  the 
particular  distributions  of  charge  and  current  in  the  beam.  Still,  one  expects  the  leading  order  loga¬ 
rithms  in  the  expressions  for  lg  and  lE,  Eqs.  (A-27,28),  to  be  correct. 

Using  now  the  zero  order  fields,  Eqs.  (3,4),  in  Eq.  (2)  we  may  write  the  condition  for  radial  force 
balance  as 


i  +  —  /a  K  -  n,o»o  +  — -S-  k  -  o 

y0 


d. 

>0  rl 


(5) 


where 


>o 


2  e 
vr£n  o — 2 
me 1 


1  “W 


(6) 


and  where  a>b  is  the  beam  plasma  frequency,  (4» rn0e2/my0)in.  Here  and  below  ftz0  retains  the 
definition  assigned  to  it  in  Eq.  (1). 

Equation  (5)  is  a  quadratic  equation  for  the  circulation  frequency,  90 ■  The  solution  which 
approaches  ft2o  as  v/yQ  —  0  is,  to  first  order  in  v/yo- 


=  ft 


zO 


1  “  ~ — V  k  +  k 

yo  or 


(7) 


where  a  -  ftz0ro/c.  Self  field  effects,  represented  by  the  v/y0  term,  are  seen  to  reduce  the  single  parti¬ 
cle  circulation  frequency  below  that  expected  for  a  zero  density;  the  correction  term  can  be  significant 
(20-30%)  in  presently  contemplated  devices.  The  general  result,  Eq.  (7),  will  be  needed  below  in  the 
derivation  of  the  first  order  equations  of  motion. 


III.  FIRST  ORDER  EQUATIONS  OF  MOTION 

In  this  section  the  equations  governing  the  motion  of  a  beam  and  motion  of  an  electron  within 
the  beam  are  obtained  and  discussed.  We  shall  consider  in  detail  only  motion  transverse  to  the  toroidal 
magnetic  field,  assuming  that  all  fields,  both  self  and  applied,  are  independent  of  9. 
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The  equations  of  motion  for  a  particle  in  the  fields  of  (A-25,  26)  to  first  order  in  the  displace¬ 
ments  from  the  reference  orbit  (r0, 0),  are  derived  in  Appendix  B.  They  are 


As  they  stand  Eqs.  (8a)  and  (8b)  are  not  easily  solved  since,  before  they  can  be  solved  for  the 
coordinates  of  a  particle  (r,,  z,)  the  beam  position  (Ar,  A z)  must  somehow  be  known  as  a  function  of 
time.  However,  a  set  of  consistent  equations  for  beam  and  particle  motion  may  be  obtained  by  per¬ 
forming  an  ensemble  average  of  Eqs.  (8a,b)  over  initial  particle  coordinates  and  velocities.  Denoting 
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such  an  average  by  brackets  it  may  be  shown  that,  as  long  as  the  beam  is  assumed  not  to  kink  (A  r,  A r 
independent  of  9),  we  will  have 

<r,>  -  A r,  <8r>  -  <8r>  -  <8r>  -  0  (10a) 

<zj>  -  A z,  <8z>  -  <8z>  -  <8z>  -  0.  (10b) 

Upon  performing  this  averaging  procedure  on  Eqs.  (8a,b)  we  will  obtain  equations  governing  the 
motion  of  the  center  of  the  beam.  These  may  subsequently  be  subtracted  from  the  original,  unaver¬ 
aged  Eqs.  (8a,b)  to  obtain  equations  governing  the  motion  of  a  single  particle  within  the  beam.  Both 
resulting  sets  of  equations  may  be  summarized  by  the  following  single  set: 

x  +  <u2jr  -  (l„0y  +  y  ft#0>>  +  F  (11a) 

y  +  u}y  -  -  Cl90x  -  y  ri8(>x  (lib) 

where  the  various  quantities  appearing  in  Eqs.  (lla,b)  are  defined  in  Table  1. 


Table  I  —  Definition  of  Quantities  Appearing  in  Equations  of  Motion,  Eqs.  (1  la,b). 


Beam  Equations 

Particle  Equations 

(x,y) 

yo/2(Ar,  Az) 

yo/2(8r,  8r) 
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nio 

_ v_ 

n 

l-n--4r 

;|^+H 
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Equations  (lla,b)  are  our  basic  starting  points  for  the  analysis  to  be  presented  below.  In  the  fol¬ 
lowing  sections  we  will  derive  and  study  the  WKB  solutions  to  Eqs.  (1  la,b).  First  we  make  a  few 
remarks  on  the  equations  themselves. 

The  term  proportional  to  x  on  the  left  hand  side  of  Eq.  (11  a)  and  the  term  proportional  to  .v  on 
the  left  hand  side  of  Eq.  (lib)  represent  radial  and  vertical  focussing  forces  respectively.  In  general  the 
coefficients  of  x  and  y  in  these  terms  are  not  equal  which  suggests  that  an  initially  circular  beam  may 
not  remain  circular.  The  value  of  n  which  makes  these  terms  equal  (the  value  required  to  maintain  a 
circular  beam  cross  section)  is 

"cir=  7  1  -  —  4  '*  +  <B  (12) 

i  yo  |  a 

which  depends  on  y0  and  therefore  on  time.  In  what  follows  we  will  leave  n  arbitrary,  though  we  shall 
assume  implicitly  that  its  value  is  close  to  «cjr.  This  is  necessary  for  self  consistency  since  we  obtained 
the  beam  self  fields  Eqs.  (A-26)  assuming  a  circular  beam  cross  section. 


In  the  case  of  constant  fields  Eqs.  (lla.b)  are  elementary.  For  this  case  we  have 


x  _  F/w2  ±  „ 
y  ~  0  +  2*.  J 


/toy-  to}  iu,l 

l - p  I 

f<uj  -  w2 


where  the  eigenfrequencies  (frequencies  of  betatron  oscillations)  are  given  by 


<0,  =  ± 


( +  u>y  +  n«o  ±  [(&>2  +  toy  +  nHy2  —  4<«) 2oj2]i/2 
—  *  1 


and  where  the  Ch  j  =  l,  2,  3,  4  are  constants. 


Stability  conditions  result  in  the  usual  way  by  requiring  u>2  >  0.  We  postpone  examination  of 
these  conditions,  however,  until  the  following  section.  We  note  here  only  that  for  values  of  y  above  a 
value  dependent  on  geometry  ( rb,a,r0,n )  but  nor  on  beam  density,  the  self  field  contributions  to  w 2  and 
<o2  fall  off  as  yo '>  rather  than  y0-3.  For  whole  beam  motion  the  value  of  y  at  which  the  r/y0  terms 

become  comparable  to  the  rft2  nja1  term  can  be  modest  (y  —  10)  for  typical  laboratory  parameters 
(rb  -  1  cm,  a  -  10  cm,  r0  -  100  cm,  n  -  0.5). 
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The  particular  solution  in  Eq.  (13)  represents  physically  for  a  particle  motion  a  first  order  radial 
shift  of  a  particle  which,  while  located  initially  at  the  reference  orbit  (r0, 0)  does  not  have  the  correct 
energy  to  be  maintained  there  by  the  local  vertical  magnetic  field.  It  therefore  moves  in  or  out  slightly 
depending  on  the  sign  of  the  energy  mismatch.  If,  however,  the  radial  focussing  forces,  represented  by 
w2,  happen  to  vanish  the  behavior  becomes  secular  (no  equilibrium  radius  exists). 

The  solution  to  the  homogeneous  part  of  Eqs.  (lla,b)  also  becomes  secular  when  a>2  =  0.  In 
fact,  when  a>2  =  0  and  <o2  ^  tux  ( n  ^  n*),  the  point  w2  =  0  corresponds  to  a  turning  point  (transition 
from  stable  to  unstable  behavior)  in  the  WK.B  solution  presented  in  the  next  section.  Since  to2  for  par¬ 
ticle  motion  will  pass  through  zero  during  acceleration,  it  becomes  important  to  examine  the  behavior  of 
the  solutions  to  Eqs.  (lla,b)  for  time  dependent  fields.  In  general,  for  slowly  time  varying  fields,  a 
numerical  solution  to  Eqs.  ( 1 1  a,b)  over  the  entire  acceleration  cycle  is  prohibitive  since  the  numerical 
integration  time  step  must  be  small  compared  to  which  in  turn  is  extremely  small  compared  to 
typical  acceleration  times.  An  explicit  solution  for  this  case  is  therefore  essential. 


IV.  MOTION  OF  BEAM  IN  SLOWLY  VARYING  EXTERNAL  FIELDS 

A.  Stability  Considerations 


If  the  coefficients  of  the  derivatives  of  x  and  y  in  Eqs.  (11  a,b)  are  slowly  varying  during  a  period 
of  a  betatron  oscillation,  the  equations  may  be  solved  by  the  WKB  method.  (See  Appendix  C.)  To 
leading  order  the  solution  is 


(x|  1  4.  A)  (<«2  —  a/2) 1/2  #/«,*’  ['  , 

w  ~  £  zr  (-)  -  «i)'" f  J  d' 


F(t')  (15) 

where  the  eigenfrequencies  are  those  given  in  (14)  in  which  now  all  quantities  may  depend  on  time. 


Kx  </./’) 

KAu') 


<ox  =  [(o>2  +  to2  +  n„20)2  -  4tojw2]1/4,  06) 

and  where  the  kernels  and  Kv(t,t')  are  given  in  Appendix  C.  The  Aj ,  j  -  1,  2,  3,  4  are  con¬ 

stants  in  this  approximation. 
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This  solution,  Eq.  06),  is  valid  far  from  any  turning  point,  i.e.  where  any  toy  vanishes.  Turning 

points  will  occur  if  to2  to2  -  0  and  if  ^  to2.  (See  below.)  Initially  we  shall  confine  attention  to  a 

cold  beam  (no  longitudinal  momentum  spread)  for  which  the  particular  solution  in  05)  vanishes  ident¬ 
ically.  Later  we  shall  comment  on  the  effect  of  temperature. 

The  solution  is  unstable  (exponentially  growing)  in  time  for  such  times  that  Im(w,)  <  0  for  any 
/  Unstable  behavior  will  occur  therefore  whenever  either  of  the  following  conditions  is  violated: 

to2to2  >  0  (17a) 

a.?  +  o,,?  +  >  2(o»2a>2),/2  (17b) 

For  n  =  „cir  (oj  2  =  a,,?)  inequality  (17a)  is  trivial  and  (17b)  gives  the  simplified  stability  condition: 

n„20  >  max(0,  -  4o>  2)  (18) 

If  n  nc„  then  both  conditions  (17a,b)  must  be  simultaneously  satisfied  for  stability.  Condition  (17a) 

in  particular  cannot  always  be  satisfied.  At  injection  ns  is  typically  quite  large  and  both  o>2  and  to,2  for 

particle  motion  (and  perhaps  for  beam  motion)  are  negative.  During  acceleration,  as  y0  increases  ns 

decreases  (ns  —  yo 3)  and  to 2  and  to2  change  sign  (for  different  values  of  y0,  if  „  ^  ,;cir);  an  instability 

"gap"  therefore  exists  while  to 2  and  to 2  have  opposite  signs. 

It  is  important  to  point  out  that  w 2  and  to2  for  beam  center  motion  (Re:  Table  I)  may  start  out 
and  remain  positive  throughout  the  injection-acceleration  cycle  while  to2  and  to2  for  particle  motion 
change  sign.  We  recall  from  Table  I  that  the  small  quantity  ( rb/a  )2  multiplies  in  the  expressions  for 
to2  and  a,,2  for  beam  center  motion  but  not  for  single  particle  motion.  Therefore  unless  is  extremely 
large  initially,  beam  center  motion  will  remain  stable. 


The  inequalities  Eq.  (17a-b)  are  illustrated  graphically  in  Fig.  3.  The  stable  regions  of  the 
2 


tOy 

2 

11  #0 

acceleration  both 


tOy 

fl#o 


and 


n 


«0 


for  particle  motion  are  negative  and  in  region  I.  In  this  region 


the  toroidal  magnetic  field  is  essential  for  stability  (modified  betatron  regime).  Following  acceleration 
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both  and  — L-  are  positive,  i.e.,  in  region  II  in  which  the  toroidal  field  is  no  longer  required 

|  J  |  J 

for  stability  (conventional  betatron  regime).  Only  by  passing  precisely  through  the  origin  (e.g.,  trajec¬ 
tory  b  in  Fig.  3)  can  instability  be  avoided  altogether.  While  the  size  of  the  instability  gap  does  not 
depend  on  the  magnitude  of  the  value  of  Im(<o7)  in  the  gap  does  and  is  inversely  proportional  to 
B„0.  Therefore  by  choosing  a  sufficiently  large  toroidal  field  it  should  be  possible  to  pass  through  the 
instability  gap  safely  (within  a  few  growth  times,  or  less). 


We  may  be  quantitative  for  a  case  in  which  toroidal  effects  may  be  neglected:  When  Eq.  (17a)  is 
violated  and  if  K?l  then  for  the  unstable  mode,  from  Eq.  (14), 


Im  <u ,  = 


yf—OliOl 


b»y 


n 


#0 


which  has  a  peak  value,  assuming  only  y0  and  not  B„0  is  changing  in  time,  of 

<d£|i„ ' 1/21  s '• 


(19) 


(20) 


If 


f  2dt  Im  ai j  =  f  JImoiy-7-^  «  1  (21) 

where  /,  and  t2  are  the  times  at  which  the  instability  gap  is  entered  and  exitted,  respectively,  then  one 
expects  that  the  transit  through  the  gap  will  not  significantly  disrupt  the  beam;  Eq.  (21)  translates  into  a 
constraint  on  y0: 


^  »  T  ft,o 4r("  -  1/2)2  (22) 

To  ■* 

If  the  acceleration  is  fast  enough  to  satisfy  Eq.  (22)  particle  motion  will  be  essentially  unaffected  by 
passage  through  the  gap.  It  should  be  possible  to  choose  a  machine  design  (i.e.,  a  sufficiently  large 
toroidal  field  and  a  field  index  close  to  1/2)  so  that  Eq.  (22)  is  well  satisfied. 


The  instability  which  occurs  while  <  0  has  an  interesting  dynamical  origin.  Let  us  consider 
the  equations  of  motion,  Eqs.  (Ila.b),  taking  F-  0,  and  taking  the  external  fields  to  be  constant  in 
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x  +  wjx  =  n  HOy  (23a) 

y  +  (o*y  =  -  ftw0x  (23b) 

These  equations  are  just  those  governing  the  motion  or  a  particle  in  an  effective  electric  field 


£'ff  =  —  io*x  (24a) 

e 

Ef  -  ~  i o},y  (24b) 

and  a  magnetic  field  B„ o/yQ.  Converting  to  polar  coordinates  p,  <t>  we  have 

reff  =  EL  p[w2  cos2  ^  +  2  sjR2  (25a) 

e 

£Jff  =  -y  p  lo>2  -  <o2l  sin</>cosif).  (25b) 


The  particle  behavior  may  be  understood  as  follows.  Let  us  assume  that  //  >  1/2,  from  which  it  fol¬ 
lows  that  a>v2  >  oil  always,  and  let  us  consider  first  the  modified  betatron  regime  (w2  <  0,  <u}.  <  0). 
£*ff  in  this  regime  is  everywhere  negative  thereby  giving  rise  to  a  clockwise  E  x  B  drift,  assuming  B„0  is 
positive.  £^ff,  which  is  much  smaller  in  magnitude  than  £*ff,  gives  a  radial  drift  of  alternating  sign  as 
the  particle  moves  from  quadrant  to  quadrant,  thereby  producing  an  elliptical  orbit.  Stable  motion  is 
established  by  balancing  the  outward  radial  electrostatic  +  outward  centrifugal  forces  against  the  f '  x  B 
confining  force. 


In  the  conventional  betatron  regime  to 2  >  0,  to2  >  0  and  the  sign  of  £*ff  is  reversed.  Azimuthal 
particle  drift  is  now  counter-clockwise  and  the  major  axis  of  the  elliptical  orbit  is  rotated  by  90°.  Stable 
motion  is  achieved  by  balancing  the  inward  radial  electrostatic  force  against  the  centrifugal  force;  the 
toroidal  field  is  no  longer  needed. 


In  the  instability  gap  £'ff  has  zeroes  at  polar  angles  given  by 


I  W  v 

COS2  0-1 - 2 


(26) 


at  which  points  the  azimuthal  drift  velocity  vanishes.  The  radial  drift  velocity,  cEln/B„,  cannot  also 
vanish  at  the  same  point.  Consequently  the  particle  drifts  radially,  with  increasing  velocity,  since 
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—  p,  at  the  angle  <£0>  as  long  as  <  0.  Increasing  the  toroidal  B  field,  thereby  reducing  the 
radial  drift  velocity,  reduces  the  growth  rate  of  this  instability,  a  fact  reflected  in  Eq.  (19). 

Typical  orbits  during  transit  of  the  instability  gap  are  illustrated  for  a  simple  case  in  Figs.  4  and  S 
in  which  results  of  a  numerical  integration  of  Eqs.  (lla,b)  are  plotted.  In  Fig.  4  condition  (21)  is  not 
well  satisfied.  The  dramatic  drift  direction  reversal  and  instability  are  evident.  In  Fig.  5  condition  (21) 
is  well  satisfied  («  is  near  1/2);  particle  motion  is  virtually  unaffected,  except  for  the  reversal  of  drift 
direction,  by  passage  through  the  gap.  The  two  graphs,  in  Figs.  4  and  5  differ  only  by  the  value  of  n 
used;  all  other  external  parameters  and  total  integration  time  are  identical. 

So  far  no  mention  has  been  made  of  the  effect  of  temperature,  the  inhomogeneous  term  in  Eqs. 
(lla,b),  on  particle  orbit  behavior  in  or  near  the  instability  gap.  Particles  having  an  energy 
mismatch— either  too  little  or  too  much  energy  to  be  maintained  at  the  reference  orbit  by  the  local 
vertical  field— will  seek  out  their  new  equilibrium  orbits  about  which  they  will  execute  betatron  oscilla¬ 
tors.  Secular  behavior  is  expected,  as  discussed  earlier,  when  vanishes. 

The  effect  of  energy  mismatch  on  a  particle  orbit  is  illustrated  in  Fig.  6  where  the  particle  of  Fig. 
5  has  been  given  an  energy  mismatch  of 

PHX  -  <P„i> 

— - —  =  y,  -  <y,>  -  0.10. 

mr0c 

The  effect  is  twofold.  The  orbit  center  shifts  slightly  outward  and  the  amplitude  of  betatron  oscillations 
following  passage  through  the  instability  gap  has  increased  by  a  factor  of  —35  over  the  zero  mismatch 
case.  Such  a  large  expansion  of  the  particle  orbits  cannot,  in  fact,  be  reliably  computed  using  the 
linearized  Eqs.  (lla,b)  used  here.  One  non-linear  effect  in  particular,  namely  the  reduction  of  beam 
density  during  the  orbit  expansion,  will  clearly  speed  the  passage  of  a  particle  through  the  instability 
gap.  (Recall  that  n,  is  proportional  to  density.)  Due  to  this  density  reduction  the  actual  degree  of  orbit 
expansion  to  be  anticipated  in  a  real  device  is  likely  to  be  significantly  less  than  that  seen  in  Fig.  6. 
Still,  these  calculations  suggest  that  a  fairly  cold  beam  will  be  required  for  successful  acceleration 
through  the  instability  gap.  Poorly  "matched*  particles  are  likely  to  be  lost  as  goes  through  zero.  It 
should  be  pointed  out  as  well  that  a  strong  toroidal  field  greatly  reduces  the  effects  of  energy  mismatch. 
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The  computer  runs  necessarily  employ  a  very  modest  toroidal  field  (600  gauss  in  the  case  of  Figs.  4-6) 
due  to  time  step  considerations.  A  stronger  field,  by  further  restricting  radial  motion,  is  expected  to 
improve  the  confinement  properties  of  a  warm  beam. 

B.  Adiabatic  Behavior 


Let  us  next  briefly  consider,  using  the  solutions  to  the  equations  of  motion,  Eq.  (15),  the  effects 
on  the  particle  orbits  of  the  removal  of  the  toroidal  magnetic  field.  The  toroidal  field  may  need  to  be 
removed  in  order  to  facilitate  beam  extraction  though  this  may  not  be  essential.  Let  us  assume  that  Eq. 
(15)  is  valid  throughout  the  acceleration  cycle,  i.e.  that  kJ  and  u,!  pass  through  zero  simultaneously 
and  that  the  solution  to  the  homogeneous  equation  (the  sum  in  Eq.  (15))  dominates  the  solution.  This 
is  certainly  true  for  matched  particles  (P„  =  <P„>  -  0)  when  it  =  1/2  and  when  toroidal  effects  may 
be  neglected  (i*/y  «  1).  One  may  show,  using  Eq.  (15)  for  such  a  case,  that  for  beam  center  motion 
in  either  the  fast  or  slow  oscillation  mode 


l(Ar)2  +  (Az)2), 

|  1  'a2 

|7  -  7"’ 

B;] it  + 

1 

l(Ar)2  +  (A.*)2], 

I  1  'A 

1 2  ■  T7"’ 

B:l  + 

f 

while  for  particle  motion  about  the  beam  center6 


(27) 


l(8r)2  +  (8z)2], 
l(8r)2  +  (8J2], 


2  “ 


B:l  +  ~B}0l 


1 

T  ~ 


B?o  +  ~B] o| 


1/2 


(28) 


where  the  subscripts  i  and  ./  correspond  to  any  initial  and  Anal  states.  The  latter  expression,  Eq.  (28), 
may  be  intc  preted  as  the  fractional  change  in  beam  cross  sectional  area.  Note  that  for  large  the 
area  of  the  orbits  ~-B„\  as  expected. 


Expressions  for  these  ratios  in  the  case  that  toroidal  effects  are  not  negligible  and  it  *  1/2  may  be 
obtained  from  Eq.  (15).  The  expressions  are  complicated,  however,  and  will  not  be  cited  here. 
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As  a  numerical  example  we  consider  a  1  kA  beam  of  1  cm  initial  radius  in  an  initial  state 
corresponding  to  y,  -  7,  B:0l  -  120  g,  Be0l  -  1.5  kg  and  a  final  state  with  yf  -  100,  BzQ/  -  1.7  kg, 
and  B„0j  -  0.  In  such  a  case  Eq.  (27)  gives  for  the  orbital  area  ratio  a  value  of  0.63  while  Eq.  (28) 
gives  for  the  ratio  of  beam  cross  sectional  areas  a  value  of  0.60. 

We  conclude  that  it  should  be  possible  both  to  accelerate  the  beam  and  to  remove  the  toroidal 
field  to  facilitate  beam  ejection  without  causing  either  the  beam  orbit  or  individual  particle  orbits  to 
expand  without  limit. 

V.  CONCLUSIONS 

The  beam  in  a  modified  betatron  can  be  stably  confined  both  during  the  acceleration  phase  and 
during  the  subsequent  gradual  removal  of  the  toroidal  magnetic  field  prior  to  beam  ejection.  As  the 
beam  is  accelerated,  however,  unless  very  special  conditions  are  satisfied,  a  region  of  instability  will  be 
passed  through;  however  if  the  time  of  transit  through  this  instability  gap  is  small  compared  to  the  time 
specified  in  Eq.  (20)  the  net  effect  should  be  small. 

As  the  toroidal  field  is  removed  to  facilitate  beam  extraction  following  acceleration  no  further  ins¬ 
tability  gaps  occur  but  the  magnitude  of  the  beam  betatron  oscillations  will  change  adiabatically.  By 
arranging  that  the  ratios,  Eqs.  (28,29),  be  near  one,  one  expects  the  beam  to  be  well  behaved  during 
the  removal  of  the  toroidal  field. 

It  should  be  remarked  however  that  changing  the  toroidal  field  changes  the  "tune”  of  the  betatron 
which,  in  general,  will  necessitate  the  passage  through  orbital  resonances  as  the  toroidal  field  is 
removed.  These  resonances,  due  to  the  periodic  encounter  by  a  particle  of  a  field  error  or  "bump"  are 
currently  under  investigation.  It  is  anticipated  that  a  condition  governing  the  minimum  speed  with 
which  B„  must  be  removed,  expressed  as  a  function  of  the  magnitude  of  the  field  error,  will  be 


obtained. 
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Appendix  A 

FIELDS  IN  THE  MODIFIED  BETATRON 


In  this  appendix  we  calculate  the  fields  seen  by  a  particle  in  a  modified  betatron.  The  particle  is 
assumed  to  be  close  to  the  axis  of  the  torus,  that  is,  the  coordinates  of  the  particle  are  taken  to  be 
(refer  to  Figs.  1  and  2  in  the  text) 

(r,z)  =  ( r0  +  A  r  +  Sr,  Az  +  Sz) 

and  all  fields  will  be  calculated  to  first  order  in  A r,  Sr,  A z,  and  Sz.  Fields  will  be  given  in  the  (r,  0,  r) 
coordinate  system  of  Fig.  1  and  all  will  be  assumed  to  be  independent  of  0.  Superscripts  a  and  s  will  be 
used  below  to  denote  applied  and  self  fields,  respectively. 

Part  I  (Applied  Fields) 


Magnetic  Field 


The  usual  weak  focussing  betatron  field  has  r  and  z  components.  The  z  component  is  taken  to 
behave  near  r0  as 

Ar  +  Sr  I 


Bf=  B:0(rJ r)n  ~  B:0  1  -  n- 


ro 


(A-l) 


where  B;0  depends  only  on  lime  and  n,  taken  as  a  constant  to  this  order,  is  the  so-called  vacuum  field 
index.  The  radial  field  is  obtained  by  requiring  (V  x  B)„  =  0  and  Br(z  -  0)  -  0  (making  the  z  -  0 
plane  a  plane  of  symmetry).  The  result  is: 


B?  =  -  nB:  c 


A z  -1-  Sz 
'0 


(A-2) 


The  applied  toroidal  field  generally  falls  off  as  r~ 1  across  the  minor  cross  section  of  the  torus: 


B»  B„ ©  1 1 


Ar  -F  Sr 
''o 
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where  B„0  depends  only  on  time.  However,  in  the  equations  of  motion  B„  multiplies  only  r  and  z  terms 
which  are  already  first  order.  Therefore  the  gradient  of  B„  does  not  enter  the  linearized  equations  of 
motion  and  we  take  only  the  zero  order  value, 

B$  =  B, o  (A-3) 

Electric  Field 

All  applied  electric  fields  are  inductive.  The  field  is  governed  by  the  changing 

central  flux  and  is  taken  to  be  a  specified  function  rtf  .'&* 

£$-£,*  (A-4) 

E„  is  negative  for  electron  acceleration  with  B:0  positive. 

Changing  the  toroidal  magnetic  field,  B„0s  will  indue  e  a  poloidal  electric  field,  the  r  and  z  com¬ 
ponents  of  which  are  easily  found: 

£?-  -  -£;B»0(Sz  +82)  (A-5) 

E?  -  -l~B"0(Ar  +  8r)  (A-6) 

2c 

where  a  dot  indicates  a  time  derivative. 

Part  II  (Self  Fields) 


Since  we  neglect  beam  diamagnetism  and  the  possibility  of  a  change  in  self  flux  due  to  time  vary¬ 
ing  beam  current  we  take  B$  =  ££  -  0.  It  remains  to  calculate  the  r  and  z  components  of  the  beam 
self  electric  and  magnetic  fields. 

Consider  a  beam  circulating  inside  a  perfectly  conducting  toroidal  chamber  of  circular  cross  sec¬ 
tion  as  shown  in  Fig.  A-l.  (The  beam  displacement  is  exaggerated  for  clarity;  we  will  assume  A  «rb). 
The  chamber  major  and  minor  radii  are  r0  and  a  respectively.  The  beam  major  and  minor  radii  are  Rb 
and  rb  respectively.  We  must  calculate  the  fields  inside  the  beam  (p  <  rb),  assuming  the  chamber  is  a 
perfect  conductor.  To  proceed  we  define  a  scalar  potential  <t>(p,4)  and  a  magnetic  flux  or  stream  func¬ 
tion  'F  (p,  <t>)  =  rA„  where  A„  is  the  usual  vector  potential.  The  equations  for  4>  and  ¥  are 
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J__d_  d* 

p  dp  p  dp 


1  d2* 
p2  dt2 


-  4 iren0(p,4>)  - 


d<t>  1  .  d<t> 

cos  4>— - sin  0— 

dp  p  d<p 


Rb  +  p  cos  t 


1  d  d'P  1  d2'!'  4ir  ,  „  ,  ,  w  /  ,  \  , 

—  —  P~  +  "TTTT  - - —  (/?<,+  p  cos t)J»(p.t)  + 


,  d'P  1.  .  d'P 

cos  4>- - sin  t  — — 

dp  p  d«£ 


(A-7) 


(A-8) 


P  dp  f  dp  )  p2  d*2  c  K  Rb+p  cos* 

where  n0,  the  beam  number  density  and  the  beam  current  density,  are  assumed  to  have  been 

specified.  Here  we  shall  take  both  n0  and  J9  constant,  independent  of  p  and  t 


The  boundary  conditions  on  <J>  and  V  are  the  same;  they  both  must  vanish  at  the  surface  of  the 
chamber,  specified  by 

p  ~  a  —  A  cos  (t  —  t).  (A-9) 

correct  to  first  order  in  A/a. 


Scalar  Potential  and  Electric  Field 

The  general  solution  for  d>,  including  the  first  toroidal  correction,  is 

P  <  rb 

p  >  rb  (A-10) 

where  q  =  -en^irr}  and<t>0-  C\  and  D'are  constants. 


d>  = 


<t»o  +  _  p2/ /’e,2)  +  cost  +  A -2-  sin  t  +  B^-  cost 

'  '  4  Rbrf  rb  rb 

<J>0  ~  2 q  In  p/rb  +  In  p/rb  cos  t  +  M  +  C  —  I  sin  t 
Kb  (  rb  P  I 


+  \b'2-  +  D  — 
'h  P 


cos  t 


Applying  now  the  correct  boundary  conditions  both  at  the  beam  surface  and  the  wall  determines 
all  of  the  constants: 


2q\na/rb 

(A-lla) 

,  A  rb  . 

A  —  -  2q - sin  0 

(A-llb) 

a  a 

B'  -  -  q-£-  lna/rb  - 

J\-2^  '»  cos* 

(A-llc) 

Kb 

ARba2  a  a 
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C'  =  o 

D'=  -2.JL. 

4 


(A-lld) 
(A  lie) 


Using  this  result  in  Eq.  (A- 10)  we  may  calculate  the  r  and  z  components  of  Es  inside  the  beam,  to  first 


£/  - - —  cos  <f>  +  —  -r—  sin  <f> 

op  p  oq> 


2q  c  .  rb  .  ,  Q  ,  a 

=  -f  8r  +  —  \r  +  -f  -  In  — 
rb  a  Rb  rb 


_s  d<t>  ■  ,  1  d4»  , 

El  =  -  —  sin  <f> - —  cos  <t> 

op  p  op 


(A-121 


=  ~  Sz  +  rAz 

r  >  /)  * 


(A-13) 


where  (8r,  8z)  =  p  (cos  <£,  sin<£)  and  (Ar,  Az)  =  A  (cos  ,  sin  «/» ) . 


Magnetic  Flux  (or  Stream)  Function  and  Magnetic  Field 


The  general  solution  for  V,  including  the  first  toroidal  correction,  is 


*0+  d'  -  P1/  rb  \  ~  t  ^~T  cos  +  A  -P-  sin  4>  +  B  cos  d>  p  < 
'  1  A  Kh  ri  rh  rh 


_  j Qjl 


I'Po  -  2 Q  In  p/rb  -  In  p/rb  cos  <t>  +  \A'-^~  +  C'-^-l  sin  <t> 


+  B'£- +  D'—  cos  <t> 
rb  P 


P  >  rb 


(A-14) 


where  Q  =  irrbJ„Rb/c  =  -  irr£enffiQRb,  /30  —  Fw(/c,  and  'Fq.  A,  fl,  A  ’,  B\  C\  and  D'  are  constants. 


Applying  the  boundary  conditions  gives 


Vo  -  20ln  — 
rb 

A  “  A  ’  —  —  2Q^-  —  sin  i li 
a  a 


(A- 15a) 


(A- 15b) 
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The  integral  over  9'  may  be  expressed  in  terms  of  the  complete  elliptic  integrals' 

4 J»  r  i 


_  rWjv'  dz,  - 

c  j«.  ..  (r„rV+u_lVw 


(1  +  -)K(-m)  -  -£(-m) 
m  m 


(A-20) 


where 


m 


4  rr' 


(r  -  rV  +  (z~  z')1  ‘ 

In  the  beam  interior  m  is  large.  Using  the  asymptotic  expansions  for  K  and  £  one  may  show  that 

(1  +  —  )K  (-m)  — — E(-m )  ~  — m_l/2(lnm  +  4ln2  -  4).  (A-21) 

mm2 

Using  Eq.  (A-21)  in  Eq.  (A-20)  the  resulting  integrals  are  elementary.  The  result,  for  the  vector 
potential  inside  the  beam  is 


At  =  - 


2  In 

»R„ 

elfb 

+  1  —  p2/  r£  -f  -  cos  $ 


—In  — — ^  +  3 

rb 


(A-22) 


where  /  -  from  which  it  follows  that  the  fields  inside  the  beam,  to  first  order  in  p,  are 


Bf 


B!=  — 


sin  <f> 


rb  |n  °Rb  _  2  P  cos  d» 
Rb  rb  rb 


(A-23) 

(A-24) 


We  may  summarize  all  of  the  foregoing  results  as  follows: 


The  applied  fields  are. 


Bf  =  -  nB.0  -Z  -- 
'0 

(A-25a) 

B»  “  B„0 

(A-25b) 

B?-B:0  jl-n 

A  r  +  8r)| 

(A-25c) 

II 

£,a  -  -  ~  B„ o  (Ar  +  8z) 

(A-25d) 

£#  “  £#o 

(A-25e) 

'Handbook  of  Mathematical  Functions,  M.  Abramowiti  and  I.  Stegun,  eds.  Dover  Publications,  ch  17. 
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£?-  S#o  (A/-  +*r) 

where  5,0,  and  E„0  are  taken  to  be  prescribed  functions  of  time. 


(A-250 


The  self  fields  are: 


Bf  =  —  2nn0efi0  |8z  +  —  Az| 

Bi  =  0 

B '!  =  2nn0ep0  Is r  +  -yA r  -  -y-  lB 
1  a  2r0 


, 2  r2 

Sr  +  -yAr  +  /c 
o2  2r0 


Ef  =  —  2nn0e 


£»  =  0 


£7  =  —  27r  «0t> 


,2 

8z  +  — y  A. 
a1 


where 


~ 


In  —  +  2  if  circulation  frequency,  0,  is  constant  across  the  beam 

rb 

In  —  +  1  if  current  density  is  constant  across  the  beam 

rb 


In  —  if  density  is  constant  across  the  beam. 
h 


(A-26a) 

(A-26b) 

(A-26c) 

(A-26d) 

(A-26e) 

(A-26f) 


(A-27) 

(A-28) 


For  times  long  compared  to  the  time  it  takes  the  magnetic  field  to  diffuse  through  the  chamber 
wall  the  result  (A-24)  shows  that  one  must  replace  a  in  the  logarithm  in  the  definition  of  lB  by  (8  r^e) 
=  2.9r0.  This  suggests  that  fields  in  an  actual  device  may  have  to  be  programmed  in  time  to  compen¬ 
sate  for  this  extra  change  (reduction)  in  B.,  in  order  to  hold  the  beam  in  place. 
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Appendix  B 

LINEARIZED  EQUATIONS  OF  MOTION  FOR  A  PARTICLE  IN  THE  MODIFIED  BETATRON 


in  this  appendix  the  equations  of  motion  for  a  particle  in  the  fields  (A-25)  and  (A-26)  are 
obtained,  correct  to  first  order  in  small  quantities. 


The  complete  equations  of  motion  are 


~(yr)  -  yr(t 2  =  — — 
dt  m 

— 

r  dt  m 


Er  H - (rtiB.  —  zB„) 

c 


Eh  +  -(ifl,  -  rB.) 


~(yz)  -CrBn  -  r»Br) 

dt  m  ■  c 


(B-l ) 

(B-2) 

(B-3) 


We  consider  first  the  linearization  of  Eq.  (B-2).  This  equation  has  an  exact  first  integral,  assum 
ing  the  fields  do  not  depend  on  W;  it  is  the  canonical  angular  momentum 


P»=r\ 


myrH - A„ 

c 


(B-4) 


We  now  write  all  quantities  Q  as  Q  -  Q0  +  Qi  where  Q|  «  and  Q0  refers  to  quantities 
evaluated  at  the  reference  orbit  (r.  r)  *  (ro,0).  Defining  •»  rii  it  is  straightforward  to  show  from 
Eq.  (B-4)  that 

V  _  !ji  -  Sl 
mydro  yd 

where 

r2 

B-a*  **  B.-o  ~  n  n<\e0o~  (g 
ro 

and  where  we  have  used  y,  -  t/MoYdf/* i/c2-  Now,  using  the  expression  for  f)0  in  Eq.  (7),  one  obtains 


«o- 


eB:\ S 


(01 


my0c 


(B-5) 
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where 


mytfro  >6  >o  a1 


(B-6) 


v 

Vo 


VO 


2  e 
— 2 
me1 


1  <*»!'* 


=  4irn0e2/my0 
a  =*  Cl:0ro/c. 


The  expression  Eq.  (B-6)  will  be  needed  next  in  the  linearization  of  the  radial  equation,  (B-l). 
Carrying  out  a  straightforward  linearization  of  Eq.  (B-l),  using  the  zero  order  fields  from  Eqs.  (A- 
25b,c)  and  Eqs.  (A-26c,d)  gives 


r  -  — 


"»y0 

+  Ki«,oyo2 


V„o 

£rl  +  —BzX 


+  n90ii  ~  —r\  ~  »of| 
Vo 


1  -  — 
Vo 


1  +  To" 


if  +  fa 


where 


(B-7) 


ft#o  “  eBJ my qc. 


Using  now  Eq.  (B-6)  and  keeping  terms  only  to  first  order  in  v/ya  and  using  Eqs.  (A-25c,d)  and 
Eqs.  (A-26c,d)  to  write 


my0 


Ve0 

Er)  +  —  Bzl 


2my0cZ 1  2y02 


Sr  + 


2H 


+  rt0on7o'’i 


we  obtain  our  final  result  for  the  radial  equation: 


(B-8) 


r\  +  ~ri  +  ft  20 
ro 

eB 


>o 


1  -  n - —  1(3  -  n)-4  fa  +  (2  -  n)  fallr, 


«o 


+  ft 


2my0c 
Pn 


*i  +  *Wi  + 


n,ft  20|8r  + 


70 


yomro 


1  -  — 
>0 


v  1 1  +  yo' 


-2 


fa  +  fa  I 


where  n,  =  a>l/2y$(lf0  is  the  "self  field  index." 


(B-9) 
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The  analogous  linearization  of  the  z  equation  (Eq.  (B-3))  is  completely  straightforward,  using  the 
fields  (A-25a,b,f)  and  (A-26a,b,f).  The  result  is 


yo 

yo' 


z  i  +  —  Z|  +  nCl 1 


zO 


1 - - - y  k  +  /fl 

Yo  a2 


2  my0c  1 


-  nMr,  +  w,n  20  8r 


(B-10) 
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Appendix  C 

WKB  SOLUTION  OF  EQUATIONS  OF  MOTION 


The  linearized  equations  of  motion  are  given  in  the  text,  Eqs.  (lla,b).  Below  we  shall  obtain 
first,  an  approximate  solution  to  the  homogeneous  version  of  Eqs.  (lla,b),  assuming  that  all 
coefficients  are  slowly  varying.  We  will  then  give  the  solution  to  the  full,  inhomogeneous  equations. 


The  homogeneous  equations  are 

x  +  m2x  -  fW  +  (C-la) 

y  +  °>}y  -  -  -  yiW  (c-ib) 

All  coefficients,  w2,  <u^,  and  Clm  will  be  assumed  to  vary  significantly  only  over  a  slow  time  scale. 
To  carry  out  a  formal  asymptotic  expansion  then  we  define 


r  =  f/X 

where  X  is  a  large  dimensionless  parameter.  Denoting  ~~  by  a  prime  ('),  Eqs.  (C-la,b)  become 

dr 


x"  +  K2(0*x  -  xn#0y'  +  yfijo’y 
y"  +  k2to 2y  -  -  xn#o*'  -  yfiTO'x. 

Now  writing 

x  -  0|(r; 

y  -  fl2(r  ;x)e'A/"Mrfr' 

we  proceed  to  express  and  a2  in  formal  asymptotic  series: 


a,(r;  X) 
a2(r;  X) 


n-0  A" 

n-0  A" 


(C-2a) 

(C*2b) 


(C-3a) 

(C-3b) 


(C-4a) 

(C-4b) 
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We  must  now  find  the  ain%  a2„,  and  <o. 

Substituting  Eqs.  (C-4a,b)  in  Eqs.  (C-2a,b)  one  finds  the  leading  order  (A2)  result 


(fc»j  -  <O2)fl|0  ~  /«ufl»oOjo  —  0 
iot +  (to*  —  <u2)a2o  ~  0 

from  which  it  follows  that  u>  must  be  one  of  the  four  quantities: 


+  v  +  ft  2 


00 


((‘“i  +  <^y  +  ft#20  ~  4<tf>2J 


1/2 


1/2 


(C-5a) 

(C-5b) 


(C-6) 


The  next  order  (A1)  relation  may  be  written,  after  some  manipulation,  as 


2(o  — 


at2  —  w  2 


a  20  +  '] 
2 


1  O)2  -  ttiy  ft  'rt) 


ft 


00 


»20 


“ft0O  +  77— (w2  -  <U2) 
“00 


o'\0  + 


I  cj'  (w2  -  to2) 

2n#°  u,  ft„0 


(C-7) 


Using  Eq.  (C-5a)  or  Eq.  (C-5b)  and  Eq.  (C-7)  an  equation  for  just  ol0  (or  just  a20)  may  be  obtained. 
The  solutions  are 


aio  “  Ai k~i/2(oj2  —  w2)l/2(w2  +  w2  +  ft«o  —  2to2)-|/2 
a2o  =  <4to~l/2(to2  —  ci>2)1/2(io 2  +  a>2  +  ftj0  —  2o)2)-l/2 
where  A  is  an  arbitrary  complex  constant. 


(C-8a) 

(C-8b) 


Using  Eq.  (C-8a,b)  and  the  definition  of  w  we  may  write  the  leading  order  WKB  solution  to  Eqs. 
(C-la,b)  as 


where  the  sum  extends  over  the  four  values  of  w  in  Eq.  (C-6)  and  where  the  Aj  are  constants. 

This  solution  is  expected  to  be  valid  as  long  as  wr,  w„,  and  ft„0  are  slowly  varying  compared  to 
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where  «  is  «ux.  u>y,  or  fl#0-  The  solution  is  therefore  expected  to  fail  when  o»y  =  0,  that  is,  near  a 
turning  point.  From  Eq.  (C-6)  this  can  happen  when 

—  0.  (C-10) 

Equation  (C-9)  is  bounded,  however,  if  in  addition  -  <uy  -  0.  Breakdown  of  Eq.  (C-9)  (and  a  tran¬ 
sition  to  unstable  behavior)  occurs  only  if  »  0  and  <n*  ^  to^. 


Once  the  solution  to  the  homogeneous  equations  have  been  found  the  solution  to  the  inhomo¬ 
geneous  equations  follows  by  the  usual  variation  of  parameters  or  some  similar  method.  Writing  four 
independent  solutions  to  the  homogeneous  equations  as 


,0)1  J 


1,  2,  3,  4 


(C- 1 1 ) 


one  finds  in  a  straightforward  way  that  a  particular  solution  to  Eqs.  (1  la,b)  in  the  text  is  given  by 


where 


(C-12) 


Kx(u')  ~-±fjktm  *0 Ht)x(k)(i')yU)U')y<m){i') 

Ky(t,t')  -  -  -jp  ejklm 

W  -  fJklmxiJ)(t)  xik)(i)y{nU)ylm)(l), 

and  where  the  summation  convention  is  understood.  The  Wronskian  If"  is  a  constant,  independent  of 
time;  its  value  is  determined  once  a  choice  is  made  for  the  xiJ\  yiJ). 
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Fig.  I  —  Cutaway  view  of  modified  betatron  geometry 
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PARTICULAR 

ELECTRON 

POSITION 


MINOR  CROSS 
OF  TORUS 


SECTION 


Coordinates  of  beam  and  particle  in  modified  betatron  Center  of  beam  is  at  lr 

at  (r.r)  "  Up  +  Ar  +  fir.  Ac  +  Sr) 


:)  -  (/„  +  A r.  Ac).  Electron  is 


30 


2 

", 

n-o 

plane.  Shaded  regions  are  stable.  Trajeclors  a  and  c  pass  through  unstable  regions. 
Only  trajectories,  such  as  b.  avoid  all  unstable  behavior 
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2 


DRIFT 

DIRECTION 


).3 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

-0.3  gr _ ^  0 


Fig.  5  —  Panicle  trajectory  (8r  vs.  fir)  in  the  modified  betatron  during  transit  of  the  instability  gap  All  parameters  are  as  in 

Fig  A  except  n  =  0.51 
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Fig.  6  —  Particle  trajectory  (8:  vs.  8r)  in  the  modifled  betatron  during  transit  of  the  instability  gap,  including  energy  mismatch. 
All  parameters  are  as  in  Fig.  5  except  an  energy  mismatch  of  (/>H|  -  <P,i>)/mr(jt  -  0.10  has  been  introduced. 
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